Determining the Energy Barrier for Decay out of Superdeformed Bands 
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An asymptotically exact quantum mechanical calculation of the matrix elements for tunneling 
through an asymmetric barrier is combined with the two-state statistical model for decay out of 
superdeformed bands to determine the energy barrier (as a function of spin) separating the su- 
perdeformed and normal-deformed wells for several nuclei in the 190 and 150 mass regions. The 
spin-dependence of the barrier leading to sudden decay out is shown to be consistent with the de- 
crease of a centrifugal barrier with decreasing angular momentum. Values of the barrier frequency 
in the two mass regions are predicted. 



I. INTRODUCTION 

Since their first experimental observation [l| , superde- 
formed (SD) nuclear states, with their strong ellipsoidal 
deformation and special set of shell closures, have offered 
a tantalizing and unique window into subatomic physics. 
Their rapid decay-out, in particular, has been the sub- 
ject of great interest (e.g^ Refs. [2-15]). In the standard 
theoretical approach [2, "SI , this process is modeled by a 
two-well potential function of deformation: Here, the nu- 
cleus is a single quantum mechanical particle, which tun- 
nels between the two wells, and can escape the system via 
electromagnetically induced decay from either. Because 
the barrier between the SD and normal-deformed (ND) 
wells is a direct consequence of nucleon-nucleon interac- 
tions, an understanding of its shape for various nuclei and 
angular momenta would be of considerable importance to 
the study of nuclear structure. Thus, a common objec- 
tive of theoretical studies is to bridge the gap between 
measured experimental data, such as lifetimes and nu- 
clear spins, and the shape of this barrier. In this Letter, 
we show that a rapid decrease in barrier height with de- 
creasing nuclear spin explains the SD decay mechanism. 

It was previously shown [10| that an elegant, two-state 
model ilij of SD decay-out is sufhcient to give an excellent 
picture of the system's time-evolution. One of the prin- 
cipal advantages of such a straightforward technique was 
the extraction from experiment of such important quanti- 
ties as the tunneling matrix element V and the spreading 
width for tunneling through the barrier . The purpose 
of the present Letter is to move beyond these phenomcno- 
logical quantities, and extract the height of the barrier 
itself as a function of nucleus and nuclear spin, as shown 
in Fig. [T] Previous approaches [lB| used a semiclassical 
treatment that did not allow for an accurate computa- 
tion of the tunneling rate prefactor. As we shall see, this 
can introduce a potential error of several orders of mag- 
nitude in the estimation of the tunneling width. In this 
paper, we go beyond previous treatments by computing 
the tunneling rate in a systematic and controlled fashion 
using a functional integral approach p7| . 
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FIG. 1: Calculated tunneling action S = 3.6W/u} versus 
spin for several SD decays in the 190 and 150 mass regions. 
W is the height of the energy barrier, and —u its curvature. 
Arrows indicate cases for which only an upper or lower bound 
on Fm has been measured. 



II. PATH-INTEGRAL APPROACH TO 
TUNNELING 

Absent additional information about the nuclear shell 
potential, the smoothest potential describing both the 
tunnel barrier and the SD well (which lies far above the 
ND yrast line at decay out) is a cubic polynomial: 



U{x) = - Xx 



(1) 



where a; is a coordinate describing the quadrupole de- 
formation of the nucleus (x = corresponds to the 
bottom of the SD well, while the ND well occurs for 
X > xb — A4uj'^/3X), uj is the oscillator frequency of 
the SD well, and M. is the inertia of the quadrupole vi- 
brational mode. With a suitable choice of the param- 
eters MllP' and A, U{x) provides a maximum entropy 
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(least-biased) fit to the unknown nuclear potential bar- 
rier. Note that for this simple potential, the barrier fre- 
quency LOB = ^/\U"{xb)\/M = LlI. 

The use of Euclidean complex-time path integrals over 
the tunneling coordinate provides a systematic approach 
for determining quantum tunneling rates at arbitrary 
temperature and dissipation fl7|, as an asymptotic ex- 
pansion in h. This method allows calculation not only 
of the leading-order exponential dependence of the tun- 
neling rate on potential parameters, but also the more 
computationally difficult subdominant asymptotics (i.e., 
the prefactor term). For tunneling through the poten- 
tial barrier ([T]) at zero temperature and friction [18-21], 
the mean-square tunneling matrix element out of the SD 
yrast state is found to be: 



36 W\ 
5 hu! J 



(2) 



where W = U{xb) = M^cu^/biX'^ is the barrier height 
(as measured from the bottom of the SD well), and Dn 
is the mean level spacing in the ND band at the energy 
of the decaying SD state. The action S to tunnel out of 
the SD state through the barrier is proportional to the 
argument of the exponential function in Eq. 
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where the factor of two is due to the power of V. 

To make contact with experiment, the tunneling ma- 
trix element may be estimated using the two-state model 
of SD decay 0, [13, [lal , which assumes the decay-out pro- 
cess is dominated by coupling of each SD state with its 
nearest-lying energy level in the ND band. The branching 
ratios and Fg = 1 — Fj^ for decay out and intraband 
decay, respectively, are determined by three rates [10] : 



F, 



N 



r^yry (Fat + r^) 
rs + rwri/(rw + ri)' 



(4) 



where Ts/h and Tn/K are the electromagnetic de- 
cay rates of the SD and ND states, respectively, and 
rV?! = / »^^-,\ is the nucleus' net tunneling rate 

through the barrier, with T = ^{Ts -t-r^v), and A the en- 
ergy difference between the SD and ND states [12] . Given 
the experimentally determined branching ratios and the 
electromagnetic widths, is known [l^: 
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The tunnelingmatrix element V may then be determined 
statistically [lB| , assuming the SD and ND levels are 
uncorrelated, and that the ND levels obey the Wigner 
surmise. The mean-square tunneling matrix element is 
found to be [la 



where a numerically negligible correction, whose relative 
size is 0{T/Dn)^, has been omitted. 

Eqs. ([2]) and ([6|) may be combined to yield an expres- 
sion for the tunneling width in terms of the properties of 
the nuclear potential barrier: 
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Note that this result for the net tunneling width, which 
includes tunneling and electromagnetic decay on an equal 
footing, differs by a factor of SF/Djv from the bare tun- 
neling width into an infinitely broadened, fully contin- 
uum ND spectrum. From the values of F5, Fjv, and Dn 
listed in Table U one sees that usage of such a bare ND- 
continuum result [2^ could result in an error of several 
orders of magnitude. 

The energy barriers (in units of the barrier frequency) 
obtained by solving Eq. ([7]) for all SD decays for which 
the four parameters, F^^ Fg, Fjy, and Dn, are known 
are listed in Table U (see also Fig. [Ij. Also listed is 
the tunneling action S = 3.6W/uj, which is a charac- 
teristic measure of the opaqueness of the barrier [1, [Isj . 
Note that W/hcu depends only weakly (logarithmically) 
on the barrier frequency lo. In the literature, the value 
huj = 0.6MeV has been used [3, [2^, but we shall deter- 
mine Lu self-consistently in Section [llll 

We note that, of the four parameters in Eq. ([7]), only 
Fs and F5 are directly measured experimentally; typi- 
cally, these are known to within a few percent. Fjv and 
Dn must be calculated theoretically, with models fit to 
experimental data. The uncertainties in F^v and Dn 
could thus be appreciable, perhaps as large as a factor 
of two or more. A better determination of these quanti- 
ties is a worthy goal for future studies of SD nuclei, but 
goes well beyond the scope of the present Letter. 

For almost all decay-out sequences, we find that the 
barrier height decreases with decreasing angular momen- 
tum. One exception is the odd-spin ^^*Hg sequence, for 
which the two highest-spin calculated barriers are so close 
that statistical fluctuations about the mean-square ma- 
trix elements of Eqs. ([2|) and ([6|) are sufficient to reverse 
the trend. This could occur, for example, due to an 
accidental near-degeneracy of the SD and ND states in 
^^'*Hg(15), which would lead to a larger than expected 
branching ratio F/v- The other exception is the first pa- 
rameter set for i9*Pb(10). The primary difference be- 
tween the first and second parameter sets for ^^^Pb is 
Ref. [l2|'s revised treatment of the pairing gap; it is thus 
seen that this consideration may play an important role 
in understanding the spin dependence of the SD decay- 
out barrier 1331. 



III. CENTRIFUGAL TUNNEL BARRIER 

Finally, we address whether the decrease in the tunnel 
barrier with decreasing spin (cf. Fig. [Hand Table|l| is con- 
sistent with the centrifugal barrier of a spinning nucleus. 
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TABLE L Barrier height W and tunneling action 5* for all SD decays for which sufficient data (branching ratios, Fs, Fjv, and 
Dff) are known. The rightmost column gives the sources of the experimental inputs as well as the sources of the estimates of 
Tn and Dn- The values of were calculated using Eq. ([5]), as discussed in Ref. [l^. The barrier frequency lj was determined 
self-consistently in the 190 and 150 mass regions, respectively, assuming the angular momentum dependence of the barrier 
height can be fit to that of a centrifugal barrier (see text). 
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^''Hg(12) 


0.26 


0.128 


0.613 
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0.049 


0.24 
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6.5 


[13, 24] 


"2Hg(10) 


0.92 
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0.37 


0.24 
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i«2pb(16) 
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0.24 
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i«2pb(14) 


0.02 
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0.24 
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i«2pb(12) 


0.34 
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0.24 
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0.24 
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"^2pb(8) 
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"'Hg(12) 


0.40 


0.086 
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0.24 
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If the superdeformed nucleus and the saddle configura- 
tion at the top of the energy barrier are described as rigid 
rotors with moments of inertia Xs f^nd Xb , respectively, 
then the barrier height W{I), as a function of the angu- 
lar momentum quantum number /, is simply the 7 = 
barrier M^(0), plus the rigid-rotor rotational increase in 
the energy of the barrier configuration, minus the rota- 
tional increase of the bottom of the SD well (from which 
W{I) is measured), i.e., 



Wil) = W{0) 



1(1 + 1) 
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Although the rigid-rotor model is a simplification, nev- 
ertheless the decrease in the barrier height between suc- 
cessive SD states can be rigorously expressed in terms of 
the kinetic moments of inertia: 



W{I) - W{I - 2) = h'^{21 - 1) 




(9) 



The kinetic moments of inertia Xg'^ of several SD yrast 
states in the 150 and 190 mass regions have been mea- 
sured. For i^^Dy^ ^ Sbh^/MeV and the aspect 
ratio 77 = b/a — 2.0 [1^, with a and b the smaller 
and larger radii of the nucleus, respectively. For ^^^Hg, 
1^ ^ 9o;i2^MeV and the aspect ratio 77 = 1.65 Q. 

The moment of inertia of the barrier configuration X^^' 
is not measured, but must be determined theoretically. 
This could be done by applying the Strutinsky shell cor- 
rection method to the cranking model (ssj . However, 
to account for pairing, we employ a phenomenological 
two-fiuid model (36j in which only the region outside the 
largest possible central sphere contributes to the moment 
of inertia. Within this two-fluid model, we find that the 
moment of inertia is 37] 



rrinTQ I — 



5/3 3 



(10) 



where the nucleus has been taken as a prolate ellipsoid of 
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revolution with aspect ratio 77 and atomic mass number 
A, rrin is the mass of a nucleon, and rg = 1.27fm. With 
these parameters, the measured kinetic moments of in- 
ertia of ^^^Dy and ^^^Hg at decay-out are reproduced to 
within 1 To leading order in the quadrupole deformation 
parameter e (see Ref. [11]), Eq. (fTO|) gives 

IW « ^A^/^m,,rl (11) 
o 

For ^^^Dy, the barrier occurs at an aspect ratio of 
r] = 1.7 [3^, so that I^g^/l's^ = 1.3 and AW = 
W{28) - W{26) = 0.21MeV. Assuming a constant bar- 
rier frequency, and comparing to the results from Ta- 
ble U {AW/hui = 0.37), implies a barrier frequency 
huj = 0.56MeV. 

For "'^^^Hg, the barrier is estimated to occur at an as- 
pect ratio of ?7 « 1.4 [s^l, so that I'^p /l'^^^ — 1.5 and 
AW = W{12) - W^(IO) = 0.12MeV. Assuming a con- 
stant barrier frequency, and comparing to the results 
from Table U {AW/fiUj = 0.51), implies a barrier fre- 
quency fito = 0.24MeV. Because the logarithmic depen- 
dence of on w almost completely cancels out in such 
a calculation, the differences AW/uj and AS" are nearly 
independent of the choice of w. 

IV. CONCLUSIONS 

In conclusion, we have determined the barrier height 
W and tunneling action S for decay-out of a superde- 
formed band by combining an asymptotically exact quan- 
tum tunneling calculation with a two-state dynamical 



model. The Table presents our numerical results for all 
superdeformed decays for which sufficient experimental 
data are known. We find that the tunnel barrier decreases 
significantly with decreasing spin during the decay-out 
process. The spin-dependence of the barrier is explained 
quantitatively in terms of the variation of the centrifu- 
gal barrier within a two-fluid model of nuclear rotation, 
which in turn allows us to self-consistently predict the 
tunnel barrier's curvature. 

The results presented in this Letter complete the 
chain of reasoning needed to connect the intriguing phe- 
nomenology of the decay-out process in superdeformed 
nuclei with an understanding of the underlying nuclear 
structure. Our results indicate that the rapidity and 
universality of the decay-out profiles can be explained 
straightforwardly within our two-state dynamical model 
by the decrease of the centrifugal barrier between the 
super-deformed and normal-deformed energy wells with 
decreasing spin. 
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